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1. Introduction 

The study of the intersections of stable and unstable manifolds of maps and flows has a strong 
influence on dynamical systems. In particular, the existence of a transverse intersection is 
associated with the onset of chaos, and gave rise to the famous horseshoe construction of 
Smale. The Poincare-Melnikov method [f36l l34l USl is a widely used technique for detecting 
such intersections. Given a system with a pair of saddles and a degenerate heteroclinic or 
saddle connection between them, the classical Melnikov function computes the rate at which 
the distance between the manifolds changes with a perturbation. 
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There have been many formulations of the Melnikov method for two-dimensional maps 
or flows [HIl [H [8l |24l and for higher-dimensional symplectic mappings [|4l|9l|2l. Recently, 
the geometric content of Melnikov 's method was exploited in order to detect heteroclinic 
intersections of Lagrangian manifolds for the case of perturbed Hamiltonian flows [|38l . Here 
it was shown that the heteroclinic orbits are in correspondence with the zeros of a geometric 
object, the so-called Melnikov one-form. 

For maps, the Melnikov function is an infinite sum whose domain is a saddle connection 
between two hyperbolic invariant sets. As usual, a simple zero of this function corresponds 
to a transverse intersection of stable and unstable manifolds of a perturbation of the original 
map. 

Melnikov's method can also be used to compute transport fluxes. In particular, a 
resonance zone for a two-dimensional mapping is a region bounded by alternating segments of 
stable and unstable manifolds that are joined at primary intersection points [|30l[T2l . Because 
the intersection points are primary, a resonance zone is bounded by a Jordan curve and has exit 
and entry sets [|T3l . The images of these sets completely define the transport properties of the 
resonance zone. Moreover, the integral of the Melnikov function between two neighbouring 
primary intersection points is the first order approximation to the geometric flux escaping 
from the resonance zone [|29l[T9l . 

The method has also been applied to the case of periodically time-dependent, volume- 
preserving flows [35] and more generally to volume-preserving maps with fixed points ll25l 
and invariant circles [|27ll . Volume-preserving maps provide perhaps the simplest, natural 
generalization of the class of area-preserving maps to higher dimensions. Moreover, they 
naturally arise in applications as the time-one Poincare map of incompressible flows — even 
when the vector field of the flow is nonautonomous. Thus the study of the dynamics of 
volume-preserving maps has application both to fluids and magnetic fields. 

Our goal in this paper is to develop, based on the theory of deformations, a general, 
geometrical description of the Melnikov displacement and to compare our theory to classical 
results. Deformation theory was first introduced in the theory of singularities [|40ll . but was 
soon used in the contexts of volume and symplectic geometry. Its application to dynamical 
systems in [|22ll2TH 7ll provide results that are close to our goals. 

Let /f be a smooth family of diffeomorphisms such that the unperturbed map /o has a 
saddle connection E between a pair of compact r -normally hyperbolic invariant manifolds. 
Let v(E) = r^M/ r S be the algebraic normal bundle of the saddle connection. We show that 
there exists a canonical C"^ section "D : S ^ called the Melnikov displacement, that 

measures the splitting of the saddle connection in first-order. We will prove that the Melnikov 
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displacement is given by the absolutely convergent series 

= Y.^fo)'^o = X(/o*)'^o (1) 

keZ keZ 

where JF^ is the vector field defined by = JF^ o . 

These sums do not converge in the tangent space T-^M, but only in the algebraic normal 
bundle v(E). The use of the algebraic normal bundle in the study of normally hyperbolic 
manifolds goes back to [fTTl . The fact that ([T]) always converges also addresses the question of 
how to deal with Melnikov 's method when the original pair of compact r -normally hyperbolic 
invariant manifolds (that are connected with the saddle connection) are not fixed with the 
perturbation. 

In addition, we will also show that the Melnikov displacement has a number of geometric 
properties. The main result in this direction is that any change of coordinates acts on the 
displacement by its puUback. This result will be used to obtain the natural action of any 
symmetries, reversing symmetries or integrals of the dynamical system on the displacement. 
Similarly, if the map preserves a symplectic or volume form, this gives additional structure to 
the displacement. For example, if is a family of exact symplectic maps and the normally 
hyperbolic invariant manifolds are fixed points, then we will show that there exists a function 
L : Z ^ M, the Melnikov potential, such that dL = i (V) co, where co is the symplectic 
two-form. This relation is reminiscent of the definition of globally Hamiltonian vector fields. 
When the normally hyperbolic invariant manifolds are not fixed points (or isolated periodic 
points), its stable and unstable manifolds are coisotropic, but not isotropic, and so the relation 
dL = i (V) CO makes no sense. 

We complete this introduction with a note on the organization of this paper. The general 
theory is developed in ^ In ^ we show how our theory reproduces the classical methods 
of Poincare and Melnikov. The study of exact symplectic maps and volume-preserving maps 
is contained in ^and ^ respectively. 

2. Melnikov displacement 

2.1. Deformation calculus 

In this subsection we present the deformation calculus for families of diffeomorphisms and 
submanifolds. We shall begin with diffeomorphisms by defining a vector field associated 
with a deformation. Next, we construct a vector field for the deformation of (immersed) 
submanifolds. Finally, we will combine these results to define and compute the Melnikov 
displacement. 
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In this paper, we consider smooth families of diffeomorphisms f^'.M^M, where M is 
an n-dimensional smooth manifold. Here, the term smooth family means that /e(<^) = /(<^, e) 
is C°° in both variables. The map /o will be called "unperturbed". 

Definition 1 (Generating Vector Field [1221). The generating vector field of a smooth family 
of diffeomorphisms f^ is the unique vector field J^f: such that 

^fe=:FeOfe. (2) 

If we regard JF^ as a nonautonomous vector field with time e, then the function = 
ft o f~^ represents its flow [|Tl Thm. 2.2.23]. Consequently, if is volume preserving, then 
jFf has zero divergence [|T1 Thm. 2.2.24], and if is (exact) symplectic, then JF^ is (globally) 
Hamiltonian. These geometric equivalences form the basis of the deformation calculus. 

Remark 1. The generating vector field J-'e was also called the perturbation vector field in [|25l ; 
however, we adopt the older terminology here. Sometimes, we will also refer to JF^ as the 
generator of the family . 

We will always use the convention that, given a family of diffeomorphisms denoted by 
italic letters , its generator is denoted by the same letter in calligraphic capitals. We have 
collected in the proposition below some well-known relations among generators, in which 
appear push-forwards and puUbacks. Recall that the puUback /* and push-forward /* of a 
diffeomorphism / : M ^ M act on a vector field A* : M — > TM as follows: 

rX = D/-1 ofXof = (Df-'X) o f 
f^X = Df o f-'X o = (DfX) o /-I . 

We note that = (/*)"! = (f-^)*. 

Next, our goal is to develop a first-order deformation calculus for invariant manifolds of 
smooth diffeomorphisms. We are mainly interested in stable and unstable invariant manifolds 
of r-normally hyperbolic manifolds for some r > 1, which unfortunately are just immersed 
submanifolds (not embedded submanifolds), and (not C°°). This gives rise to a few 
technicalities. Recall that a map g : N ^ M is an immersion when its differential has 
maximal rank everywhere. If g is one-to-one onto its image, then W = g(N) is an immersed 
submanifold of the same dimension as A'^. We will denote an immersion by = g(N) M 
or simply W M when the immersion g does not matter. For brevity, we will sometimes 
omit the term "immersed". 

Example 1. If V7 is the stable (resp., unstable) invariant manifold of a fixed point of a 
diffeomorphism, then = W (resp., = M"), where s and u are the number of stable 
and unstable directions at the hyperbolic point. When the fixed point is hyperbolic, n = s-\-u. 
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Stable and unstable manifolds are typically not embedded because they can have points of 
accumulation. In this case the immersion is not a proper map. 

We consider families of submanifolds of the form We = ge(N) M, where g(^, e) is 
in both variables, for some r > 1 . All the elements of such a family are diffeomorphic 
(as immersed submanifolds), because they are diffeomorphic to the same "base" manifold N. 
Just as for /o, the unperturbed submanifold is denoted by Wq. 

Definition 2 (Adapted Deformation). IfW^ M is a C family of immersed submanifolds, 
a family of diffeomorphisms (pe : Wq ^ is an adapted deformation when cpo = Idwo and 
e) is C in both variables. 

Adapted deformations exist since it suffices to take (j)^ = ge ° Sq^ ■ While there is 
quite a bit of freedom in the choice of (pe, only its normal component is relevant, since this 
measures the actual motion of We with e. The normal component will be defined using the 
algebraic normal bundle. For an immersed submanifold W M, this is defined as the set of 
equivalence classes 

v{W) = TwM/TW . (3) 

When M is Riemannian, this normal bundle is isomorphic to the more familiar geometric 
normal bundle, T (cf. [[T6l p.96]). In general v is a manifold of dimension n and is defined 
independently of any inner product structure onT M. 

Definition 3 (Displacement Vector Field). The displacement vector field of a C family of 
immersed submanifolds We M is the C^~^ section 

d 



de 



(j^eiO + T^WQ (4) 
e=0 



where (j)e '■ Wq ^ We is any adapted deformation. 

The displacement vector field is well-defined; that is, its definition is independent of the 
choice of the adapted deformation as is shown in the following lemma. 

Lemma 1. Let (j)e,4'e '■ W^o ^ be two adapted deformations. Then 

e T^Wo , V<f G Wo . 



d ~ d 
de de 



e=0 

-ii 



Proof. For each fixed <f e V^o^ the map e ^ c(e) = (p^ {cjjeii)) describes a C curve in Wq 
such that c(0) = <f and c'(0) e T^Wq. Using the fact that 4>e{0 = (peic{e)), we then have 



' d ~ d 
de de 



= DMOc'(0)=c'(0)eT^Wo, 

Je=0 

because D0o(^) = Wr^ivo- ^ 



d 

0e o /o = o /o + D/o — 

6=0 56 



•re 

6=0 
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When the submanifold is invariant under a diffeomorphism, its deformations are related 
by means of a fundamental iterative relationship between the generating vector field of the 
family of diffeomorphisms and the displacement vector field of the family of submanifolds. 

Proposition 2. Let fe be a smooth family of diffeomorphisms, and We M be a family 
of immersed submanifolds that are invariant under f^. Then 

f*V -V = f*J^o (5) 

on the unperturbed submanifold Wo, where V is the displacement vector field (|?]). 

Proof. The tangent space TWq is invariant under the puUback /J", so the term f^V is well- 
defined as a section of the normal bundle v{Wq). If (p^ '■ Wq ^ is any adapted 
deformation, then = fe ° (j^e ° fo^ is as well. Differentiating 4'e ° h = fe <Pe with 
respect to e yields 

d 

de 

where we used ([2]). However, by lemma [T] the displacement (|4]) is independent of the adapted 
deformation, so 

P o /o = J-o o /o + DfoV . 

Applying D/q~^ o /q to both sides finishes the proof. □ 
The identity (|5]) is equivalent toV = (fo)*!) + !Fq. Thus, we can work either with push- 
forwards or puUbacks. To obtain the Melnikov displacement we will iterate these identities 
on the stable and unstable manifolds of a family of diffeomorphisms. 

2.2. Normally hyperbolic invariant manifolds and saddle connections 

The Melnikov displacement will be defined for a saddle connection between a pair of normally 
hyperbolic invariant manifolds. In this section we recall the definitions of these objects. There 
are many slightly different definitions of normally hyperbolic manifolds, see [fTTl . In this 
paper, we adopt the following. 

Definition 4 (Normally Hyperbolic Invariant Manifold). Let A c M be a submanifold 
invariant under a smooth diffeomorphism f : M ^ M . We say that A is r -normally 
hyperbolic when there exist a Riemann structure on TM, a constant 1 e (0, 1), and a 
continuous invariant splitting 

TaM = El®El®TaA, Va G A (6) 
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such that if L^a^ : E^'" ■ Tf(a)A are the canonical restrictions of 

the linear map D/(a) : TaM Tf(^a)M associated to the splitting (|^, then 

(i) \^a\ < ^' and 

(ii) |L^|'|(L^)-i| <X 

for all I = 0, 1, . . . , r, and for all a e A. 

As is usual in the literature, the term normally hyperbolic will be taken to mean 1- 
normally hyperbolic. Note that setting / = in the previous definition implies that the 
linearizations and L" of / restricted to the stable and unstable spaces of A have the uniform 
bounds 



I ' 



< A < 1 , 



Va e A . 



(7) 



In this paper we will assume that each r-normally hyperbolic invariant manifold A is 
compact, although it would be sufficient that our diffeomorphisms be uniformly in some 
neighbourhood of A. We will also assume, without loss of generality, that A is connected. 
One consequence is that the sets 



W\A) = W\A, f) 



and 



= W"(A) = W"(A, /) = 



^ &M : lim dist(/^(<^), A) = 



f e M : lim dist(/^(f), A) = 

A:— > — oo 



are C immersed submanifolds of M that are tangent at A to TA^E^'^, see [fTTll . In particular, 
TaW^ n TaW^ = TA. Moreover, A and its stable and unstable invariant manifolds are 
persistent: given any smooth family of diffeomorphisms fe such that / = /o, then for each 
small enough e there exists a nearby r-normally hyperbolic invariant manifold with 
families of immersed submanifolds Wp^ = V7^'"(Af , /f). 

To compute the Melnikov displacement, we will need to show that certain series are 
geometrically convergent; the following lemma is a key component in this proof. 

Lemma 3. Let f : M ^ M be a diffeomorphism with a compact normally hyperbolic 
invariant manifold A, and fix a point ^ ^ = W^{A, f). Then given any splitting 
T^M = T^W^ ® A^^, there exists a constant ji e (0, 1) and an integer nQ > Q such that 

where : N^^ ^f„+i are the restrictions o/D/(<f„) : T^„M T^^+iM to the subspaces 
N^^ = Df"i^)[N^], = f"(.0- A similar bound holds for W". 
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Proof. There exists a unique point a ^ A such that lim„^+oo dist(^„, a„) = 0, where 
a„ = f^ia). According to the A-Lemma for normally hyperbolic manifolds [|5l, the 
complementary subspaces N^^ tend to the unstable subspaces , as n ^ +00. Thus, the 
maps : A^^ ^4+1 tend to the unstable restrictions L"^ : E^^ ^a„+i as n ^ +00, 
and the lemma follows from Q. It suffices to take any ju g (A, I). □ 
The Melnikov displacement will be defined as a function on the normal bundle of a 
saddle connection, which is defined as follows. 

Definition 5 (Saddle Connection). Let f : M ^ M be a diffeomorphism with a pair of 
compact normally hyperbolic invariant manifolds A and B. A saddle connection between A 
and B is an invariant submanifold S c W"(A) fl W^(B) such that 

Z = T^WiA) = T^W'iB) 

for all e S. 

Remark 2. The coincidence of the tangent spaces is needed in order that the manifolds have 
the same algebraic normal bundles on the saddle connection: 

v{i:) = viW\A))h=viW\B))h 
since T^M/T^X = T^M/T^WiA) = T^M/T^W\B), for all f g S. 

By definition, dim Z = dim V7"(A) = dim W^(B), and the manifolds A and B are not 
part of the saddle connection. The simplest (and most common) saddle connections are of 
the form S = W"(A) \ A = W^(B) \ 5. In this case, we say that the unperturbed invariant 
manifolds are completely doubled. Many Melnikov problems studied in the literature fall into 
this category. Nevertheless, in some problems there may exist points f e W"(A) fl W^(B), 
^ ^ AU B, such that T^WiA) # T^WiB), see [9]. In that case the saddle connection 
is strictly contained in the intersection of the stable and unstable invariant manifolds: S ^ 
(W"(A) n W'iB)) \ (A U B). 



2.3. Displacement vector fields of stable and unstable invariant manifolds 

In this subsection we use the fundamental iterative equation (|5]) to obtain infinite series for the 
displacements of the stable and unstable manifolds. These series are absolutely convergent, 
but only, as we must stress, when they are evaluated on their corresponding normal bundles. 
Indeed, the tangential components of these series can be unbounded. Consequently, in order 
to compute these sums, each term must be projected onto the normal bundle. An example will 



be given in §2.5 
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The proof of the following proposition is inspired by a proof given in [|2l, the main 
difference is that our setting is more geometric. 

Proposition 4. Let be a smooth family of diffeomorphisms such that the unperturbed 
map fo has a compact normally hyperbolic invariant manifold Aq with stable and unstable 
invariant manifolds Wq^. Then the displacement vector fields T>^ : Wq v(Wq) and 
"D" : Wq v{Wq) of the families of perturbed stable and unstable invariant manifolds 
are given by the absolutely convergent series 

k>l k<0 

Proof We prove the claim about the stable displacement V^; the unstable result is obtained 
analogously. Repeatedly applying the iterative formula ([5]), yields 

n 
k=l 

for any integer n > 1. Therefore, it suffices to check that the term (/q*)"!^^ tends 
geometrically to zero on the normal bundle of the unperturbed stable manifold. For any 
point ^ e = W'iAo,fo), let 4 ^ f^\^), and ^ V\in) ^ T^„M. Then 
{f*)"V^{^) = {LnT^Vl where L„ = T)f^{^) : T^M T^^M. We must show that 
(Ln)~^'D^ tends geometrically to zero as an element of the quotient space T^M/T^Wq. 

Given any splitting T^M = T^Wq ® N^, let : T^M N^^ be the projections onto 
the linear subspaces A^^ = DfQ(^)[N^], and let : N^^ ^i„+i be the corresponding 
restrictions of the linear maps L„ = D/o(<f„) : Tf„M 7f„+iM. Then, 



n+l 

because L„ = L„_i • ■ ■ LiLq and o = H^^^ o . Let no be the integer referred to 
in lemma [ij and define /„o = H^^q^ |(L^)~^ |. Lemma [s] gives the bound 
li-rA'/f ^-l7-)S <- / ,,"-"0 |rr^r)M < / 1 7:)M 

for some ju e (0, 1). Moreover, the sequence {T>l)n>o is bounded due to the compactness 
of the normally hyperbolic manifold Aq and the continuity of the displacement vector field 
P^(<f). Therefore, |n^ {Ln)~^'D^\ tends geometrically to zero as n ^ 0, and so, (L„)~^P^ 
tends geometrically to zero in the quotient space T^M/T^Wq. □ 
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2.4. Melnikov displacement 

We now will use the displacement vector fields to study the splitting of a saddle connection 
upon perturbation. As usual, /<= denotes a smooth family of diffeomorphisms such that the 
unperturbed map fo has a saddle connection S c M^"(Ao, /o) H W^(Bo, fo) between a pair of 
compact r -normally hyperbolic invariant manifolds Aq and Bq. These manifolds persist and 
remain r -normally hyperbolic for small e. 

We want to study the distance between the perturbed manifolds = V7"(Af , /f) and 
= W^(B(;, /f) The growth rate of this distance with e is obtained simply by taking the 
difference between the displacement vector fields of both families. 

Definition 6 (Melnikov Displacement). Under the previous assumptions, the Melnikov 
displacement is the canonical C''~^ section of the normal bundle defined by 

V = V'' -V^ : H ^ vCL) (8) 

where : Wq v(W^q) andV^ : Wq v{Wq) are the displacement vector fields ofW^ 
and W^, respectively. 

Remark 3. The Melnikov displacement V makes sense only on the saddle connection, where 
the tangent spaces of the unperturbed invariant manifolds Wq and Wq coincide. Away from 
S , the difference V^—V^ is undefined because each term is a section of a different (algebraic) 
normal bundle, see remark |2j 

Corollary 5. Let fe : M ^ M be a smooth family of diffeomorphisms verifying the 
assumptions of definition^ Then its Melnikov displacement V : S vCL) is given by 
the absolutely convergent sums 

= X(/o)'-^o = Z(/o*)'-^o . (9) 

Proof From proposition Q we get that P = - = Y^ke^ifof^Q = Z^ez(/o*)^-^o, 
where the last equality follows from the identity (/o*)^ = (/q*)"^- □ 
The Melnikov displacement has been defined in a canonical way as a section of the 
normal bundle, as such it has strong geometric properties. We will show next that any change 
of variables acts as a puUback on it. This will imply a number of geometrical properties, 
for example that the displacement is invariant under the puUback and the push-forward of 
the unperturbed map. In addition, we will see that if there exist symmetries, reversors, 
first integrals, or saddle connections, then these have natural implications on the Melnikov 
displacement. These claims are the subject of proposition [6} 
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We recall that a diffeomorphism f : M ^ M is symmetric when there exists a 
diffeomorphism s . M ^ M such that f o s = s o f, and then s is called a symmetry 
of the map /. Analogously, / is reversible when there exists a diffeomorphism r : M ^ M 
such that f or = ro f~^, and then r is called a reversor of the map /. In many applications, r 
is an involution, r^ = Id, though this need not be the case [fTSll . Finally, a function / : M ^ M 
is a first integral of / when I o f = I . 

Proposition 6. Let fc be a smooth family of diffeomorphisms verifying the assumptions of 
definition^and let V be its Melnikov displacement. 

(i) Given any smooth family of changes of variables h^ : M ^ M, the family ff^ = 
hj^ o fg o hf also verifies the assumptions of definition^ its Melnikov displacement 
T> is defined on the saddle connection £ = /jq ^(Z), and 

V = h^V . (10) 

(ii) The Melnikov displacement is invariant by the pullback and the push-forward of the 
unperturbed map. That is, 

f*V = (fo).V = V . (11) 

(Hi) If ff has a smooth family of: 

(a) Symmetries s^ '■ M ^ M such that sqCL) = S, then 

s*V = (so).V = V . (12) 

(b) Reversors r^ : M ^ M such that roCL) = Z, then 

r*V = {ro).V = -V . (13) 

If in addition, ro is an involution and its fixed set Rq = e M : ro(<f) = is a 
submanifold that intersects S transversely at some point ^o, then T>{<^o) = 0. 

(c) First integrals 1^ : M ^ R such that U a I~^({0}), then 

i(V)dIo^dIo{V) = 0. (14) 

(d) Saddle connections c fl such that Sq = £, then V = 0. 

(iv) If there exists a vector field X : M ^ T M such that itsfiow commutes with f^, then the 
Lie derivative of the Melnikov displacement with respect to X vanishes, that is, 

LxV = . (15) 

It is very important to stress that the above results make sense as identities on the normal 
bundle of saddle connections. The relation ( 10) makes sense because the pullback /zq maps 
T S onto r S, whereas (11) makes sense because T S is invariant by the pullback /q* and the 



push-forward (/o)*. Similar arguments apply to (12) and (13). The identities (14) and (15) 
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12 



make sense because T Z is contained in the kernel of the one-form d/o and the vector field 
X is tangent to Z, respectively. The hypothesis U 5^ c I~^({0}) means that and 
are contained in the same level set of the first integrals, which can be assumed to be the zero 
level without loss of generality. This holds, for instance, in the transitive homoclinic case: 
5f = Af and Af is transitive. 

Proof. We can write a geometric proof based on the definition of the Melnikov displacement 
as a canonical section of the normal bundle, or a computational proof using the formulae (|9]). 
We follow the geometric approach. 

(i) The first claims are obvious. For the last one, note that if : Wq^ V7|'" are 
any adapted deformations (for the family ), then = hj^ o o Hq are as well (for the 



family f^). Differentiating /if c 
HQoho + Dho 



de 



o ho with respect to e yields 

d 



6=0 



0^'" o h 







6=0 



where we used that i^/z^ = TY^ o /i^ . Hence, Hq oho + D/zqX'^'" = "P^'" o Hq, so the difference 



56 



V = — V^ verifies the relation V o ho = DhoV, which is equivalent to ( 10). 
(ii) If we take h^ = in (i), then = f^^of^of^ = f^ and Z 



fo\^) = 2- 



Therefore, V = V and (11) follows from ( 10). 



(iii.a) If we take h^ = Se in (i), then f^ = s^^ o f^ o Se = fe and 1. = Sq ^Z) = Z. 



Therefore, V = V and ( |12[ ) follows from ([lOj). 

(iii.b) If we take h^ = r^ in (i), then = r~^ o f^ o r^ 



/-I and Z 



'(Z) 



Z. Moreover, a stable (respectively, unstable) manifold of a map becomes an unstable 



(respectively, stable) for the inverse map. Therefore, T> = —V and ( 13 ) follows from ( 10). 

Next, we assume that ro is an involution whose fixed set Ro intersects Z transversely 
at (fo- That is, r^ = Id, ro(<fo) = ^o^ and T^qM = T^^Z ® T^qRq- We want to prove that 
P((fo) = in the normal bundle v (Z), or equivalently that P(<fo) ^ Tfo ^ • 

Since ro(fo) = <fo and r^ = Id, the square of the linear endomorphism Dro(<fo) : T^qM ^ 
T^qM is the identity map. This implies that Dro(<fo) is diagonalizable and its spectrum is 
contained in the set {—1, 1}, so r^yA/ = ® E~, where = ker(Dro(<fo) tW). We claim 
that £■+ = T^qRo and E~ = Tf^Z. The first claim follows from the fact that involutions 
are locally conjugate to their linear parts at fixed points. Since T^^Z is invariant under 
Z)ro(<fo) and complementary to T^qRq = in T^qM, we get the second claim. Finally, 
the evaluation of ([13]) at the point ^ = ^o, yields Z)ro((fo)^^(fo) = -^^(^oC'fo)) = -^^(6), 
and so Vi^o) ^ E~ = T^yZ. 



Canonical Melnikov theory for diffeomorphisms 13 
(iii.c) Let I, = Iq + eh + 0(e^). Obviously, Wf^' C I~\{0}). Differentiating 



o 01'" = with respect to e yields d/oCP"'") = -h. Thus, D = D" - verifies ( 14). 



(iii.d) In this case, we can take 0| = 0", and soV = — vanishes on S = Z 



0- 



(iv) Let(p' : M ^ M be the flow of the vector field: = X o (p' and 0° = Id. Thus, 
given any t, the (constant) family = (p' is a smooth family of symmetries of fe- Next, 
using relation ( 12), we get that (cp'yv = V for any t, and ( 15 ) follows by definition of the 



Lie derivative. □ 
These results have been extensively used in the literature. The invariance of Melnikov 
objects under the unperturbed map gives rise to periodicities when suitable coordinates are 
used; examples can be found in [[T4l [8ll20l. Item (iii.d) implies a simple splitting criterion: if 
the Melnikov displacement does not vanish identically, the separatrix splits [HI [371. Upper 
bounds on the number of uniform first integrals of the family can be deduced from 
item (iii.c), see [[32l[4TI . This result has also been used in to establish necessary and sufficient 
conditions for uniform integrability of analytic, exact-symplectic maps fl21il . Symmetries have 
also been extensively used, for example to improve the lower bound obtained by Morse theory 
for the number of critical points of some Melnikov potentials [0 (We will discuss Melnikov 



potentials in §4.2 ). Relation (15) is similar to Noether's theorem, since the existence of a 
continuous symmetry — the flow of the vector field X — implies a conservation law for the 
Melnikov displacement. Finally, fixed sets of reversors can be used to guarantee the existence 
of heteroclinic points and zeros of Melnikov functions; this is an old trick, see [[TOl . 

In the classical Melnikov method, one uses simple zeros of a Melnikov function to predict 
the transverse intersection of the invariant manifolds. We next show that this result also holds 
for the Melnikov displacement. 

Theorem 7. If i^q is a simple zero of the Melnikov displacement then the perturbed 
invariant manifolds and intersect transversely at some point <ff = + 0{e) for 
small enough e. 

Proof. By definition, the saddle connection S is an /o-invariant submanifold of M. Let 
n : v(E) Ebe the projection of the normal bundle onto S. There is a tubular 
neighbourhood N of 1, that is diffeomorphic to v(S) via a diffeomorphism y/ : ^(1,), 
as illustrated in figure [T] Since is an open neighbourhood of S in M, each deformation of 
Z occurs inside A^, for e small enough. 

We note that under the diffeomorphism y/, and for e small enough, deformations of Z 
can be thought as deformations of the zero section, 0;r, in v(S). Indeed, a deformation of S 
can be parametrized by a section E v(E). Let Uq c S be an open set that contains fo- 
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Figure 1. A tubular neighbourhood that is diffeomorphic to the normal bundle ;r : v ( E ) — > S . 
Let and f/" be deformations of Uq such that 

Let A| = I// (i7|) and = y/ {U^) be the images in v(S) of two deformations of the saddle 
connection, corresponding to the images of the stable and unstable manifolds. We want to 
show that A| and intersect transversely for e small enough when V has a simple zero. 

We parametrize each manifold with a section Uq — > That is, if e is small, then 

^1'" : Uo^ v(S) given by 

#"(0 = 7r-i(^)nA^'" 

are sections of the normal bundle restricted to Uq. Notice that the functions : Uq ^ N C. 
M given by 

j.s,u , —1 Is,u 
(Pe = W o (Pe 

are adapted deformations of Uq with images [/|'". In fact, both are perturbations of the zero 
section. From the definition of displacements we have that 

where 

and ), are the displacements of ^l'"^. Therefore, each manifold A| and A" is the image 

of ^1 and ^" and they intersect transversely at ^|(<fo) = ^"('fo) if and only if ^q is a simple 
zero of the section 0" — 0| . 



Canonical Melnikov theory for diffeomorphisms 15 

Now we use the standard "blow up" argument so that the implicit function theorem can 
be applied. Let : Z v(S) be the section given by 



l^l'"" for 6 # 
P^'" for 6 = 



Notice that 4''" is if (^|'" is C^ and moreover that, when f 7^ 0, is a simple zero of 
^^-4>l if and only if it is for0^-0|. Finally, since D^//(OP((f) = (^^ - (^^) (0 and and 
A" are the images of ^| and the implicit function theorem implies that if is a simple 
zero of V, then A| and A" intersect transversely near ^q, for e small enough. Thus and 
intersect transversely near ^q, for e small enough. □ 

2.5. Example: perturbed Suris map 

As a simple example, we consider the generalized standard map 

/:TxM^TxM, f{x, y) = {x + y - V'{x), y - V'{x)) (16) 

where V : T ^ T is a periodic potential. It is easy to see that / preserves area and orientation 
and that its fixed points have the form {x* , 0) where V'{x*) = 0. Such fixed points are saddles 
if and only if V"{x*) < 0, because tr(D/(;c*, 0)) = 2 - V"{x*). 

Following McMillan [|3T1l . we can find a generalized standard map with a saddle 
connection between two saddle points if we choose a diffeomorphism c : M ^ M such 
that c(jc + 2) = c~H-« + 1) + 1 = c{x) + 2, and let 

V'{x) = 2x-c{x)-c~\x) . (17) 

To see this, first note that, with this choice, the force V'{x) is periodic with period one. 
Moreover, if x* is a hyperbolic fixed point of c (that h, c{x*) = x* and < c'{x*) ^ 1), then 

V"{x*) = 2 - c'{x*) - {c-^)'{x*) = 2 - c'{x*) - \/c'{x*) < . 

Thus, (.x;*, 0) is a saddle fixed point. 

Moreover, the graphs of the functions /± (x) = x — ^x) are invariant and the dynamics 
on these sets is very simple: 

f\x, x-{x)) = {c'{x), x-{c\x))) , 

(18) 

f\x,x+{x)) = {c-\x),x+{c-\x))) 
for all ^ e Z. These graphs contain saddle connections if we choose a pair of neighbouring 
fixed points a and b of c, so that c has no fixed points in (a,b). Suppose further that a 
(respectively, b) is a stable (respectively, unstable) fixed point of c, so that limyt^+00 c^ix) = 
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a and liniyt^+oo c~^(x) = b for all x e (a,b). Then A = (a, 0) and B = (b, 0) are saddle 
points of the map /, and 

I^± = [(x,x±(x))eTxR:xe(a,b)} 

are saddle connections between them: c W(A) n and Z+ c W"(A) n 

see figure |2] 




Figure 2. Saddle connections and E_ for the Suris map with potential ( [20) i and the 
diffeomorphisms x \—> c^^(x) for = 0.2. Also shown are images of a point on the saddle 
connection obtained by iterating x t-^ c(x). 



It is known that a generalized standard map with a potential of the form ( 17 1 is typically 
nonintegrable Il27l . An integrable example, /o, for each ju e (0, 1) is obtained when the 
diffeomorphism c is given by 

2 



c(x) = Ca(x) = — arctan 

71 



an + l)tan(7rx/2) + (/i - 1) \ 
V(/i - l)tan(7rx/2) + (/i + 1)/ 
for —1 < jc < 1. This function, when extended to M using c{x + 2) = c{x) + 2, gives the 
period-one potential ( [TT] ) 

Voix) = - arctan -—^—^ df, S= ^ . (20) 

n Jo \l + d cos (In t) / (1+//)^ 

A first integral of the map is I{x,y) = cos ny + Scosn(2x — y) [|39ll33ll26ll . 

The diffeomorphism ( [T9] ) is conjugate to a Mobius transformation, and it is easy to find 

an explicit formula for its iterations: 

{c^f(x) = c^,(x), VfceZ. (21) 
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The points a = — 1/2 and b = 1/2 are fixed points of c. The point a is stable and b is unstable 
because c'{a) = e (0, 1) and c'(b) = 1/ ju > I. 



We now perturb the integrable map by modifying the potential (20): 

Ve(x) = Vo(x) + eU(x) . (22) 

We compute the Melnikov displacement using Q. For this calculation we must first compute 
the vector field J-'q, which for the generalized standard map (16) with potential (22) is 

J^o(x,y) = -U'ix-y)il,lf . 

Using the relations ( [TS] ) and (21 ), the series (|9]) for the displacement V is easily computed. 
However, in order to ensure convergence of the series, we must take into account the fact 
that it is only the normal component of the displacement that is desired; indeed, the iteration 
of the tangential component is not bounded. Let 11^ : T-^M be the canonical 

projection onto the normal bundle (|3]). Due to the fact that Z is invariant under /o, we have 
that o /q* = /q* o n^. To avoid numerical errors, we project each term in the sum (j9|. We 
display in figure |3] the projection 11^ (V) as a function of x along Z_, for two perturbative 
potentials U{x). 




Figure 3. Normal component of V for the Suris map on S_ with jU — 0.2. The solid curve is 
for U{x) — — cos(2;rx)/2;r and the dashed curve is for U (x) — sin(2;rx)/2;r. 
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The Melnikov displacement ([8]) generalizes the classical methods used to detect the splitting 
of separatrices due to Poincare and Melnikov. 

3.1. Poincare method 



Poincare's method p6l is based on the existence of first integrals, such as an energy function, 
for the unperturbed system. Assume that /o has a saddle connection S between a pair of 
hyperbolic fixed points a and b. For the simplest case, S has codimension one, and a single 
first integral suffices: /o/o = /. Saddle connections with a higher codimension can be treated 
in a similar way if there are sufficiently many integrals. The splitting between the stable and 
unstable manifolds is measured by the rate of change in the first integral / with e on the saddle 
connection. That is, we define the Poincare function M/ : S ^ M by 

d 



Ml 

de 



(locP^-IocPl) (23) 
where 0" and 0| are deformations adapted to the perturbed invariant manifolds. 



Proposition 8. Under the above assumptions, 

Ml = i iV)dI = dliV) 

where T> is the Melnikov displacement @. 

The proof is a simple computation. The term d/ (V) makes sense because the first integral 
is constant on the saddle connection, and so the tangent space T S in contained in the kernel 
of d/. If d/| 2 is nondegenerate, then a point <^o ^ 2 is a simple zero of M/ if and only if it is 
a simple zero of V. In this case, the perturbed invariant manifolds intersect transversely near 
<fo for e small enough. Therefore, the Melnikov displacement V generalizes the estimates of 
splitting in the Poincare style. 

3.2. Melnikov method 

The classical Melnikov method [34] is based on estimating the movement of a manifold in 
a direction normal to the separatrix. To define the normal, the phase space M is assumed to 
have a Riemannian inner product (•,•). As before, assume for simplicity that S is a saddle 
connection of codimension one. The appropriate normal to Z is called an adapted normal 
vector field. 
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Definition 7 (Adapted Normal Vector Field rt25l ). Let 1. be a fo-invariant submanifold of 
codimension one. A vector field tj : "L ^ T-^M is adapted normal when 

(i) rj is nondegenerate: r}{^) ^ 0, for all e Z, 

(ii) rj is normal: {r](^),v) = O,for all e E and v e 7^2; and 

(Hi) T] is invariant: /o*('/, 3^> = {f], fQy), for all vector field y : ^ Ty,M. 
Using this notation, the classical Melnikov function : Z ^ M is 

{n,re-(l^l) (24) 

6=0 

where 0" and 0| are deformations adapted to the perturbed invariant manifolds. Consequently, 
Mfj is related to the Melnikov displacement ([8]) by 

M„ = {rj,V). (25) 



d 



Since rj is nondegenerate, Myj and V have the same simple zeros. We also note that (25 1 makes 
sense because rj is normal, and so the component of V in the tangent space T S does not play 
any role. 

Remark 4. Adapted normal vector fields can be obtained as gradients of nondegenerate first 
integrals. Let / be a smooth first integral of /o such that its gradient V/ does not vanish 
on S. Given a Riemannian structure, the gradient is the unique vector field such that 
/ (3;) dl = (V/, 3^) for any vector field 3^ on M. Therefore, /q ( V/, 3^) = ( V(/ o /o), f*y), 
and f] = V/ is an adapted normal vector field. Obviously, the Poincare function ( [23] ) and the 
Melnikov function ([24]) coincide when rj = S/ 1 . 



4. Exact symplectic maps 

4.1. Basic results 

In this section, we will see how the Melnikov displacement ([8]) generalizes previous theories 
developed for hyperbolic fixed points of exact symplectic maps. We will also show why 
general normally hyperbolic invariant manifolds can not be studied in the same way. 

A 2n -dimensional manifold M is exact symplectic when it admits a nondegenerate two- 
form CO such that co = — d/l for some Liouville one-form k. The typical example of an exact 
symplectic manifold is provided by a cotangent bundle M = T* Q, together with the canonical 
forms coq = dx Ady and Aq = ydx, in cotangent coordinates (x, y). 

A map f : M ^ M is exact symplectic H §y A = §f{^y ) ^ for ^i^Y closed path y o M 
or, equivalently, if there exists a generating (or primitive) function S : M ^ M. such that 
f*A — A = dS. In particular, an exact symplectic map is symplectic: f*co = co. 
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A submanifold AT of M is exact isotropic if A = ^^^^-^ A for any closed path y (Z N 
or, equivalently, if there exists a generating function L : N ^ M. such that j^X = dL. Here, 
: N "-^ M denotes the natural inclusion map. In particular, an exact isotropic submanifold 
is isotropic: j^co = 0. The maximal dimension of an isotropic submanifold is n, and when 
the dimension is n, the submanifolds is called Lagrangian. 

A vector field T : M ^ TM is globally Hamiltonian if there exists a Hamiltonian 
function : M — > R such that 

/ (JT) (y = dH . 

We stress two key properties of the symplectic case. First, the generator of a family 
of exact symplectic maps is globally Hamiltonian and there exists a simple relation between 
this Hamiltonian and the generating function of the maps. Second, the stable and unstable 
invariant manifolds of a connected normally hyperbolic invariant submanifold A of an (exact) 
symplectic map are (exact) isotropic if and only if A is a fixed point, in which case they 
are Lagrangian (that is, n -dimensional). The second property is an obstruction to develop 
a symplectic version of our canonical Melnikov theory for general normally hyperbolic 
invariant submanifolds, we shall do it just for fixed points. 

These properties are well-known, but we prove both for completeness. 

Proposition 9. Let be a family of exact symplectic maps with generating function and 
generating vector field T^. Then, is globally Hamiltonian with Hamiltonian 

dS 

H, = - o /-I . (26) 

de 

Proof. By definition, = J^g o fg, so the Lie derivative with respect to J^g is Lj-^X = 
(/;)-i^(/;A). Using Cartan's formula: L^X = i(X)dX + di(X)X, and taking the 
derivative with respect to e of the relation dS^ = f*X — X, we get 

= f /;a = A = /; (/ (j^,) dx + di {t^) x) . 

de de 
Rearranging this yields 

/ {T,) dX = (/-I) d^ - di {T,) X = d\^-^o /-I - i {T,) A j . 

Finally, / {T^) X = X(J^f), since A is a one-form. □ 

Proposition 10. Let f : M ^ M be a diffeomorphism with a connected, normally 
hyperbolic, invariant manifold A. If f is (exact) symplectic, the stable and unstable invariant 
manifolds W^'^ = VK^'"(A, /) are (exact) isotropic if and only if A is a hyperbolic fixed point. 
Moreover, in this case the submanifolds W^'^ are (exact) Lagrangian. 
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Proof. Since A is normally hyperbolic, ([6]) implies that for each a g A, s = dimEl, 
u = dimE" and c = dimr^A sum to dimr^M = dimM = In. Since / is symplectic 
s = u = n — c/2. Therefore dim = dim = u + c = n+ c/2. Suppose that W^'^ 
are isotopic; then since isotopic manifolds have maximum dimension n, we must have c = 0. 
Consequently, dim A = and since A is connected it is a hyperbolic fixed point. 

Conversely, assume that dim A = 0, so that A = {a} for some hyperbolic fixed point a 
and dim W^'^ = n. To prove that is Lagrangian, take any two vectors u,v e T^W^. We 
know that Df^u and D/'^o tend to zero as ^ +oo, since the stable directions are uniformly 
contracted. Since / preserves co, we have 

co{u,v) = {f'^Ycoiu, v) = co{Df''u, d/o) — ^ co(0, 0) = 

ask ^ +00. 

If, in addition, / is exact symplectic, then for every closed loop y , A = §jr(^y^ ^. Now 
suppose that y c V7^ so that f^(y ) ^ {a} a.s k ^ +oo. Then 



Finally, the (exact) Lagrangian character of follows from the fact that it is the stable 
manifold for / ~ ^ □ 

4.2. Melnikov method for hyperbolic fixed points of exact symplectic maps 

Let /f : M ^ M be a family of exact symplectic maps such that /o has an exact Lagrangian 
saddle connection A between two hyperbolic fixed points a and b. Note that in this case 
that the stable and unstable manifolds are as smooth as the map /f. We assume that 



Ho(a) = Hoib), where //^ is the Hamiltonian (26). Without loss of generality, we can assume 
that Ho(a) = Ho(b) = 0. 

The natural measure of splitting for this case is a real valued function L, the Melnikov 
potential, whose derivative measures the splitting [9|. In other words, nondegenerate critical 
points of the Melnikov potential predict transverse splitting. In this subsection, we shall find 
the relation between the one form dL introduced by [[381 and the Melnikov displacement. 

Definition 8 (Melnikov Potential). For a saddle connection A and Melnikov displacement 
V : A ^ v(A), the function L : A ^ M implicitly defined by 

dL = jX{i(V)co). (27) 

is the Melnikov potential^ 

t For simplicity, henceforth we will write this relation as / (V) oj — dL. 
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We will show next that L is indeed defined by ( [27] ) and that its critical points correspond 
to zeros of the displacement. 

Proposition 11. The pullback of the one-form i (V) co to the saddle connection A is well- 
defined and exact. In particular, there exists a function L : A ^ M, determined uniquely up 



to additive constants that obeys ( [27| ). 

Moreover, the set of simple zeros of the Melnikov displacement T> coincides with the set 
of nondegenerate critical points of the function L. 

Remark 5. Let be a submanifold of M and X : N ^ v{N) a section of its normal bundle. 
Then the pullback of the one-form i {X) co to the submanifold is well-defined if and only if 
N is isotropic. This has to do with the fact that {X) co) is well-defined if and only if 

co(X -\-y,Z) = co{X,Z) 

for any vector fields X : N ^ T^M andy, Z : N ^ TN. 

Proof. The pullback is well-defined because A is Lagrangian, see remark [5j With regard to 
the exactness, it suffices to prove that <f^ jj^(i (V) co) = for any closed path y c A. Let y a 
closed path contained in the saddle connection, then 

/ jXii (V) co)=l jXii if^V) if*co)) = I {V) CO) 

Jy Jy -J foiy) 

where, since /o is symplectic, the Melnikov displacement is invariant under the pullback of 
/o, and the saddle connection is invariant under /q. Finally, since y e W^Q), fo), we obtain 
that 

/ jX(i (V) co)=<f (V) co)-^ <f r^(i (V) co) = 

Jy Jfl^(y) Jb 

as ^ -|-oo. 

The equivalence between simple zeros of V and nondegenerate critical points of L 
follows from the Lagrangian character of the saddle connection A. □ 
An expUcit series for L can be obtained using the Hamiltonian ( [26] ). 



Corollary 12. The Melnikov potential ( [27p is given by the absolutely convergent series 

L = Y,Hoof^. (28) 

Proof. Since fo is symplectic, and the generator JFq is globally Hamiltonian with Hamiltonian 
Hq: i (JFq) CO = dHQ,we deduce that 

/ {V)co = Y, i {(fof^o) CO = X(/o (-^o) co = Y, d(^o o /^) = dL . 

keZ keZ keZ 
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The series converges absolutely because HQ{a) = Ho(b) = 0, liniyt^-oo /q (f) = d, and 
limyt^+oo /o (<f) = i>. In fact, it converges at a geometric rate. □ 
While the Melnikov potential has been used many times for exact symplectic, twist, and 
Hamiltonian maps, the formulation give here, using the power of deformation theory, is more 
elegant. 

The Melnikov potential introduced here is identical to the one defined in [9|. This can be 



checked by direct comparison of the formula (28) with the formula (2.7) of the cited paper. 



using (26 1 to express Hq in terms of the derivative of the generating function 5f at 6 = 0. 

We also note that it is impossible to define a "Melnikov potential" on the saddle 
connection of normally hyperbolic invariant manifolds with nonzero dimension, because then 



the saddle connection is not isotropic (proposition 10) and so the identity (27 1 makes no sense 
(remark[5]). Nevertheless, even in this case the Melnikov displacement is defined. 



4.3. Area-preserving maps 

In this subsection, we restrict to the two-dimensional case in order to show in a simple way 
that the Melnikov potential and the classical Melnikov function are transparently related. 
Moreover, we will see that in some cases there exist geometric obstructions to the non- 
vanishing of Melnikov functions. These points are most easily seen by choosing a special 
time-like parametrization of the saddle connection. 

We consider the standard symplectic structure on the plane, (M, co) = (M?, dx A dy), 
and let J be the standard 2x2 symplectic matrix: co(u,v) = Jv = {u, Jv). 

Suppose that : ^ is a family of diffeomorphisms preserving area and 
orientation, and is the generating Hamiltonian for JF^ , the generator for fg . We assume that 
the unperturbed map has a saddle connection S c W^(a) fl W^(b) between two hyperbolic 
fixed points a and b such that Ho{a) = Ho(b) = 0. Note that the unperturbed map not need 
be integrable. 

The key point is that, on a one-dimensional saddle connection, there is a parametrization 
a : M ^ S such that 

/o(a(f)) = a{t + I) , lim a(t) = a , lim a(t) = b . 

t—> — oo f— >+oo 

Remark 6. In many cases, such parametrizations can be expressed in terms of elementary 
functions, and the Melnikov function can be explicitly computed f[T4l [8ll23ll. 

Consequently, a provides a diffeomorphism between the saddle connection S and the 
real line, so that objects defined over S can be considered as depending on the real variable t. 
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Thus, for example, the Melnikov potential ( [28] ) can be replaced by L o a to become a function 
L : M ^ M given by 



L(0 = ^Hoiait + k)) . 



(29) 



Here we abuse the notation by not giving the function a new name. 



Our goal is to show that the classical Melnikov function (25 1, or rather the composition 
M o a, can be computed by differentiating L: 

M{t) = L'it) . (30) 

In addition we will show that these functions have the properties 

• Periodicity: L{t + 1) = L{t) and M{t + 1) = M{t). 

• In each fundamental domain [t, t + 1), M must vanish. Indeed, j'^^ M(s)ds = for 
any f e M. (This property will be generalized to volume-preserving maps in the next 
section.) 

• Near each simple zero of M(t) or nondegenerate critical point of L(t) there is transverse 
intersection of the stable and unstable manifolds. 



Given (29) and (30), the first two properties are obvious. Periodicity is simply a 



consequence of the invariance of (29) under t ^ t + \ ; this invariance under the unperturbed 
map is a property of all of the Melnikov objects introduced so far. The second property is a 



simple consequence of periodicity and integration of (30). 



To show that M(t) is actually the classical Melnikov function (25 ) we must construct an 
adapted normal vector field rj on A, recall definition |7j We claim that, when thought of as a 



function of t, such a vector field rj : 
rjit) = Ja'(t) . 



is given by 



This claim is proved in lemma 13 at the end of this subsection. 



The relation between M and the classical Melnikov function ( [25| ) follows from a 
straightforward computation of the derivative using = —I: 

L'it) = ^(V//o(«(? + k)), a'it + k)) 



keZ 



^(Ja'it + k), -JVHoiait + k))) 

keZ 

Y,{riit + k),J^o(a(t+k))) . 



keZ 
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Thus 



keZ 



which is the obvious form of (25 1 under the parametrization a . This verifies (30) and the final 
property. 

To end this subsection, it remains to prove the claim about the vector field rj. 

Lemma 13. The vector field rj : 1, ^ M? defined by T](a(t)) = Ja'(t) is an adapted normal 
vector field on the saddle connection S. 

Proof. Since fo is symplectic, Df^ JDfo = J. Consequently, 

Ja'it + 1) = Df-\ait + l)fjDf-\a(t + l))a'(t + 1) 
= Df-\ait + l)fja'it). 

Thus, given any vector field y : S ^ M^, we have 

f*{rj, y)(a(t)) = {rj(a(t + 1)), y(a(t + 1))> 
= {Ja'{t + l),yia{t + m 
= {Ja'it), Df-Hait + l))yiait + 1))) 
= {Ja'it), Df-'ifoiait))yifoiait))) 
= {r},foy)(ait)). 

Moreover, since 7 is antisymmetric, {rj, a') = 0. Therefore, according to definition [v} 
rj : Z ^ is an adapted normal vector field. □ 



5. Volume-preserving maps 

For the case of a volume-preserving mapping with a codimension-one saddle connection, an 
adapted normal field formulation of the Melnikov function also applies [12511 . Here we show 
how to relate this to the Melnikov displacement ([8]). 

Let /e : M ^ M be a family of volume-preserving diffeomorphisms on an oriented 
n-dimensional manifold M with volume form Q. such that fo has a codimension-one saddle 
connection Z between two normally hyperbolic invariant sets A and B. 

We start with a simple lemma about the generator for . 

Proposition 14. Let J-'e be the generator of a volume-preserving smooth family . Then 

(i) The divergence of J-'^ with respect to Q. is zero. 

(ii) The one-form i iJ-'^) Q is closed. 
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(Hi) If M is simply connected, then i (JF^) Q is exact. 

Proof. From [1, Thm. 2.2.21], we have that Ljr Q, = 0. Since the divergence is defined by 
div(jFf)Q = Ljr Q,, this implies that the divergence vanishes. Moreover, since dQ. = and 
L^^ Q = di (J^e) Q + i (J^e) dQ, then di (J^e) Q = 0, implying (ii) and (iii). □ 
In order to find an invariant nondegenerate (n — l)-form from any adapted normal field, 
we assume that M has a Riemannian metric (•, •). 

Proposition 15. Ifrj is an adapted vector field (cf. definition^, then 

iin)^ ,,,, 

{tj, rj) 

is a nondegenerate (n — I) -form on S that is invariant under the restriction /olz- 

Proof. By definition rj is nonzero, so that co^ is nondegenerate. To prove that fQCO,^ = co^ on 
the saddle connection S , we introduce the vector field : S ^ T S defined by 

2 _ fpV fl_ 

This vector field is tangent to S, because {rj, Z^) =0. Now, we compute the difference 

_^ ^ i (/o» /o ^ i in) Q ^ i (/o» ^ i in) ^ CI 
Jo^r, CO, ^^^^^^^ ^^^^^ ^^^^^ ly ,j 

Hence, it suffices to see that the {n — l)-form i (^Z,) Q vanishes identically on the tangent 
space T S . This follows from the fact that Z, is tangent to S and dim E = n — 1 . □ 
The Melnikov function associated with rj is defined using corj. 

Definition 9 (Volume-Preserving Melnikov Function). Let V be the Melnikov displace- 
ment ([§]). Given an adapted vector field tj on S, we define the Melnikov function M, : S ^ M 
as the unique C''~^ function such that 

Mr,co, = i (V) a (32) 

as (n — \)-forms on Z. 

The previous definition first appeared in [|25ll . Note that / (V) Q. is an {n — l)-form, but 
it is possible that it might be degenerate; if this were the case then zeros of Mrj need not 
correspond to those of V. We will show next that this is not the case. 

Proposition 16. The Melnikov function is invariant under the map fo and M^ = {t},T>). 
Moreover, a point <^o is a zero ofM, if and only if it is a zero ofV. 
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Proof. Using (31 1, we obtain 

i (V) a - {t], V)co^ = i (o) Q 

where v = V rj. Since v e TZ, we conclude that / (o) Q = 0, as an (n — 1)- 

{}], rj) 

form on Z, and thus Mrj = {f},V). Since is an adapted field, and V is invariant, 
f*{fj,V) = {tj, f*V) = {tj,V). Therefore, f*M^ = M^. Finally, M,,(<fo) = only when 
V(^q) = 0, since rj is nonzero and normal to T E , and V is not in the tangent space. □ 



We would like to show, as we did for the area-preserving case in section 4.3 , that the 
volume-preserving Melnikov function necessarily has zeros on E . To do this, we will show 
that the integral of M,j with respect to the measure co^ is zero. This can be accomplished 
by dividing the saddle connection into pieces — fundamental domains — that are mapped into 
each other by /q. 

Definition 10 (Proper Boundary). Let Abe a compact normally hyperbolic invariant manifold 
of a diffeomorphism fo with stable manifold VF^(A). A proper boundary, y , is a submanifold 
ofW^{A) that bounds an isolating neighbourhood of A in W^{A). In other words, y is proper 
if there is an closed submanifold (A) such that 

(i) BW^A) = y,and 

(ii) /o(W;(A))c/nf(W^;(A)). 

We refer to the closed set Wy (A) as the stable manifold starting at y . 

Similarly, for the unstable manifold, a submanifold a c W"(A) is proper it is proper for 
/q"^ However, in this case we define the unstable manifold up to a, denoted by W^{A), as 
the interior of the local manifold that corresponds to f^^. 

Notice that the definition is not symmetric, because (A) is a closed subset of W^(A), 
while W^(A) is open in M^"(A). The asymmetry is just a technicality in order to simplify 
some proofs. 

Definition 11 (Fundamental Domain). Let A be a hyperbolic invariant set for fo. A 
submanifold with boundary P is a fundamental domain of W^{A) if there exists some proper 
boundary y e VF**(A) such that 

P = W;{A)\W%,y{A). 
Equivalently, a fundamental domain in VF"(A) is a manifold with boundary of the form 

P = W'^{A)\W''^_,^^{A), 
where a is a proper boundary in W^(A). 
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In each case, the fundamental domain is a manifold with boundary dV =y U foiy). 
An immediate consequence of the definition is that all the forward and backward iterations 
of a fundamental domain are also fundamental. It is easy to see that proper boundaries 
always exist, and in fact, the unstable manifold can be decomposed as the disjoint union 
of fundamental domains: 

W'(A)\A=\JfSiP). 

The importance of fundamental domains is that much of the information about the entire 
manifold can be found by looking only at these submanifolds. For example, as discussed in 
[|27ll . the topology of the intersections of and can be studied by restricting to P. 

Proposition 17. Let M be a simply connected manifold and f(::M^M a family of volume- 
preserving maps such that fo has a saddle connection Z with fundamental domain P. Ifco^ 
is the (n — I) -form defined in (31 ) and M,j is the Melnikov function {32), then Jp M^co^ = 0. 

Proof. The fundamental domain P is a submanifold with boundary, such that dP = y U/o(7 ), 
where y is a closed curve that does not intersect /o(7 ). If we give an orientation [P] to P , the 
induced orientation on the boundary satisfies [y ] = — [/o(7 )]. According to proposition [141 
the form i (JFq) Q. is exact. Thus there exists an {n — 2)-form H such that dH = i (JFq) Q. 
From (|9]) we conclude that 

/ {V) Q = i (Y,(fo)'^o)n = X(/o)'(' (-^o) = dX(/o = d¥ 

\keZ / keZ keZ 

where = ^^g^C/o*)^^- Finally, from definition|9| the Stokes's theorem, and the invariance 
= ¥, we get that the integral 

/ M,co,= [ i(V)Q = [ d^= [ ^ = - [ f*^ 

Jp Jp Jp JdP Jy Jy 

vanishes. □ 
The previous result implies that the stable and unstable manifolds of a perturbed saddle 
connection necessarily intersect. Examples of such intersections were computed for the case 



M = M — where the hypothesis of proposition 17 are satisfied — in 112511271 



6. Conclusion and future research 

We have studied a general theory of the Melnikov method that can be applied to many different 
settings. The formula Q for the Melnikov displacement V generalizes many of the classical 
methods that use a normal vector field to measure the displacement with respect to a natural 



Canonical Melnikov theory for diffeomorphisms 29 

direction. For example, when the saddle connection is defined as the level set of a first 
integral /, the classical Poincare function Mj = dI(V) measures the splitting as the rate 
of change of the first integral. If there is a Riemannian structure and an associated adapted 
normal vector field rj, then the classical Melnikov function = {rj,V) measures the rate of 
change of the splitting in this normal direction. For the case of exact symplectic maps with 
saddle connections between hyperbolic fixed points, the Melnikov potential L is defined on a 
Lagrangian submanifold and acts as a generator for the displacement: dL = i (V) co. 

This Melnikov theory can be extended to other situations and can be applied in many 
problems. For instance, one can study billiard dynamics inside a perturbed ellipsoid, 
following a program initiated in [|6l|3l . It turns out that the billiard map inside an ellipsoid is an 
exact symplectic diffeomorphism defined on the cotangent bundle of the ellipsoid, which has 
a two-dimensional normally hyperbolic invariant manifold with a three-dimensional saddle 
connection. Therefore, ellipsoidal billiards represent a strong motivation for a more detailed 
study of general normally hyperbolic invariant manifolds in a symplectic framework. This is 
a work in progress. 

In discrete volume geometry there are many open Melnikov problems, since the 
application of Melnikov methods to volume -preserving maps began just a few years ago [|25l 
l27l . We plan to continue this program in several ways. As a first step, we plan to obtain 
bounds on the number of primary heteroclinic orbits in terms of the degree of the polynomial 
perturbation [28] . 
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